Highly anisotropic dissipative hydrodynamics 



Michael Strickland 

Physics Department, Gettysburg College, Gettysburg, PA 17325 United States 

Abstract. The quark gluon plasma generated in ultrarelativistic heavy ion collisions may possess 
sizable momentum-space anisotropics that cause the longitudinal and transverse pressures in the 
local rest frame to be significantly different. We review recent attempts to derive a dynamical 
framework that can reliably describe systems that possess a high degree of momentum-space 
{N) anisotropy. The dynamical framework that has been developed can describe the evolution of the 

quark gluon plasma ranging from the longitudinal free-streaming limit to the ideal hydrodynamical 
O limit. 
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The Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) are 
studying the behavior of nuclear matter at high energy densities, e ^> 1 GeV/fm 3 , using 
relativistic heavy ion collisions. The goal of these experiments is not only to generate a 

^ deconfined quark-gluon plasma (QGP), but to also study its properties such as transport 

properties, color opacity, etc. One complicating factor is that the QGP generated in such 
collisions lasts for only a few fm/c and during this time the bulk properties of the system, 
e.g. energy density and pressure, can change rapidly. In addition, it is possible that in 
certain regimes the system is far from equilibrium/isotropy. Dynamical models that can 

qq reliably describe the evolution of the system on the fm/c timescale are necessary in order 

to make sound phenomenological predictions. 

One of the key outstanding questions in dynamical models of the QGP is to what 
extent is the QGP isotropic in the local rest frame. Early indications from ideal hydro- 

. £h dynamical fits to experimental data for the elliptic flow of hadrons indicated that the 

data were consistent with isotropization at x ~ 0.5 fm/c after the initial nuclear impact 
[1,2]. This rather short time scale prompted a plethora of papers attempting to solve the 
early isotropization/thermalization puzzle. In the intervening ten years ideal hydrody- 
namics has been replaced by viscous hydrodynamics [3-6] as the method of choice for 
simulating the bulk dynamics of the QGP. 

In viscous hydrodynamics the energy-momentum tensor in the local rest frame is 
naturally anisotropic with ^ where and are the local rest frame 
transverse and longitudinal pressures, respectively. The relative amount of momentum- 
space anisotropy is encoded in the shear tensor and in viscous hydrodynamics one has 
the freedom to choose, in addition to initial temperature and flow profiles, an initial value 
for the shear profile. In the last years it has emerged that phenomenological predictions 
are rather insensitive to the assumed initial shear profile, see e.g. [8]. Further studies have 



shown that the period of large momentum- space anisotropy can persist for a few fm/c 
[9]. In addition, it can be shown that momentum- space anisotropics will be particularly 
large at the longitudinal and transverse edges of the plasma where the system is dilute 
and poorly approximated by naive viscous hydrodynamical treatments [10]. 

In order to properly address the question of the evolution of systems that possess large 
momentum-space anisotropics, one needs to go beyond traditional viscous hydrodynam- 
ical treatments. Viscous hydrodynamics relies on an implicit assumption that the shear 
correction is small and that one can linearize around the isotropic ideal background. If 
the shear correction is large, a new framework is needed. This has motivated the de- 
velopment of reorganizations of viscous hydrodynamics in which one incorporates the 
possibility of large momentum-space anisotropics into the leading order of the approxi- 
mation [9, 1 1-16]. The framework developed has been dubbed anisotropic hydrodynam- 
ics. This method is capable of describing ideal hydrodynamics through free streaming 
in a single framework in which the expansion is organized around the smallness of the 
off-diagonal components of the energy-momentum tensor. 



FORMALISM 

If we have a system that is azimuthally symmetric in local rest frame (LRF) momenta, 
then the energy momentum tensor can be expressed in terms a timelike four- vector u^ 1 
and a spacelike four- vector z^ which are mutually orthogonal [15, 16] 

T^ v = {£+ ^ T )u^u v - <?> T g^ + {&> L - ^> T )z^z v , (1) 

where <S is the energy density, is the transverse pressure, and S^l is the longitudinal 
pressure. The spacelike four- vector z^ is directed along the beamline direction of a heavy 
ion collision and is the four-velocity of the LRF. 

One can derive dynamical equations for the energy density, pressures, and by 
taking moments of the Boltzmann equation p^d^f(x,p) = —C[f] where / is the one- 
particle distribution function and C is the collision kernel. The moments are defined by 
multiplying the left and right hand sides of the Boltzmann equation by various powers of 
the four-momentum and then integrating over momentum space. This can be achieved 
via the n th moment integral operator J^„ = / d% p^p^ 2 ■ ■ ■ p^ n where n > is an integer 
and 

The resulting zeroth moment of the Boltzmann equation can be compactly written as 

Dn + n0=J o , (3) 

where D = u^d^, 9 = n is the number density, and 7o = (2ft)~ 3 / <i 3 p P^f/ Po is 

a particle number source. This equation governs the evolution of the number density. In 
number conserving theories Jq = 0; however, in number non-conserving theories such 
asQCD,7 ^0. 

The first moment of the Boltzmann equation is equivalent to the requirement of energy 
and momentum conservation 

d^ v = 0, (4) 



Taking projections parallel and tranverse to one obtains 



D£ + + ± )e + (^ L -^)u v D LZ v = 0, (5) 

and 

{S + ^ L )Du a -V a ^ L +z a D L {^ L -^ 1 _)+z a {^L-^L)BL 

+(^ L - ^)D L z a - {@> L - ^ ± )u a u v D LZ v = , (6) 

respectively, where Di = z^du and 9i = z^dn. 

In order to proceed one can use symmetries to restrict the form of the one-particle 
distribution function. In the case that the system is azimuthally symmetric in LRF 
momenta and consists of massless particles, it suffices to introduce a single scale, A, and 
a dimensionless anisotropy parameter ^ [16]. These parameters define a deformation of 
an arbitrary isotropic LRF distribution function originally introduced in Ref. [17] 

/(f , x, p) = / iso ((p 2 T + [1 + £ (f , x)]p 2 L )/A 2 (t, x)) , (7) 

where we have explicitly indicated that A and t, are functions of space and time. With 
this form for the one-particle LRF distribution function and the general equations for 
the components of the energy momentum tensor obtained by taking moments of the 
Boltzmann equation, one can derive dynamical equations for A and £, using 

d 3 p n iso (A) 

U — 7z > W 



(2nr V /TTT' 

where n iso is the isotropic (t, = 0) number density. One can also evaluate the energy- 
momentum tensor in the LRF 



(27t) 3 po 



and obtain [18] 



<f(A,<*) = =<%(!;) £[ S0 (A), (10) 

^r(A^) = ^(T xx + Tyy) = M T (t;)& iso (A), (11) 

^l(A^) = -r/ = ^ L (£)^iso(A), (12) 

where ^i so and S^o are the isotropic {£, = 0) pressure and energy density, respectively, 
and 

*«,-4(l±«^-«)). 

m) = I ( « + 'ff»-' ) . (15, 




The equation of state can be imposed as a relationship between & lS0 and ^i so . 

In the case that the system is boost invariant along the beamline direction and one 
assumes an ideal equation of state, the form (7) results in 

^Dl; -6D(logA) -29 = IT (l -&/*(S)y/lTt) , (16) 

from the zeroth moment of the Boltzmann equation, where we have used the relaxation 
time approximation for the collision kernel with a relaxation rate T = 2^ ,1 / 4 (<^)A/5f] 
with rj = rj/S with rj being the shear viscosity and S being the entropy density [12]. The 
first moment in the boost-invariant case becomes 



—u 



$[3a(S)+a ± (S)]D(jj?) =0 L (Z)D 1 S+Aig 1 (S)D 1 QDgA), (17) 

u 2 

from the first moment, where = d r uo + V± ■ u±, D = uod r + -^-u± • Vj_, Dj_ = 
z ■ (u_l x Vy) = u x d y - u y d x , u_l = (u x , Uy), and = 1 + u\. 

RESULTS AND OUTLOOK 

Having obtained Eqs. (16) and (17) one can solve them to find the spatio-temporal 
evolution of the energy momentum tensor of an anisotropic system. By construction 
the transverse and longitudinal pressures obtained in the evolution are guaranteed to be 
positive. This should be contrasted to solutions to second order viscous hydrodynamics 



which can result in negative pressures at early times and near the edges of the matter 
[10, 12]. In Fig. 1 we plot the ratio of the LRF longitudinal and transverse pressures at 
three different proper times. For this figure we assumed a non-central b = 7 fm collision 
and used a Monte-Carlo sampled Glauber wounded-nucleon profile, a central b = 
isotropic temperature of Aq = Tq = 0.6 GeV at To = 0.25 fm/c, £o = 0, and Anr\ /S = 1. 

As can be seen from Fig. 1, even though the system starts out being perfectly isotropic 
at To, after a very short amount of time the system develops a pressure anisotropy on 
the order of ~ 0.4 - 0.5 that only slowly relaxes back towards isotropy. One 

should note, importantly, that this figure is generated assuming the best case scenario for 
isotropization, namely that Anr\jS = 1. If one chooses larger values of 7], then one finds 
larger momentum- space anisotropics at all times. However, regardless of the value of rj, 
when evolved using Eqs. (16) and (17), the pressures (and in particular the longitudinal 
pressure) remain positive at all times in the entire transverse plane. 

In closing I have briefly reviewed the derivation of (2+l)-dimensional anisotropic 
hydrodynamics. The full detailed derivation of the equations along with the numerical 
algorithms necessary are contained in Ref. [16]. The next step in the development of 
anisotropic hydrodynamics is to relax the assumption of azimuthal symmetry of the LRF 
one-particle distribution function in momentum-space. This represents work in progress. 
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